Weak-Coupling Approach to Hole-Doped S = 1 Haldane Systems 
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As a weak-coupling analogue of hole-doped S = 1 Haldane systems, we study two models for 
coupled chains via Hund coupling; coupled Hubbard chains, and a Hubbard chain coupled with an 
S = 1/2 Heisenberg chain. The fixed point properties of these models are investigated by using 
bosonization and renormalization group methods. The effect of randomness on these fixed points is 
also studied. It is found that the presence of the disorder parameter inherent in the Haldane state 
in the former model suppresses the Anderson localization for weak randomness, and stabilizes the 
Tomonaga-Luttinger liquid state with the spin gap. 



The Heisenberg spin chain with integer spin-shows a re- 
markable feature, so-called the Haldane gapEl. Recently, 
hole-doping into S — 1 Haldane gap systems was realized 
in Y2BaNiO50. It inspires theoretical interest in the ef- 
fects of carrier-doping into Haldane systems, which have 
not been studied well so farfl. In this paper we consider 
weak-coupling models for hole-doped Haldane systems 
(HDHS); coupled Hubbard chains via Hund coupling (re- 
ferred to as model I) and a Hubbard chain coupled with 
an S = 1/2 Heisenberg chain (model II). The model I cor- 
responds to the case where the energy levels of electrons 
which compose the 5=1 state are almost degenerate, 
whereas the model II to the case where these levels are 
largely separated, and the lower level can be considered 
to be localized. Coupled chain problems have been stud- 
ied extensively in several contextsuila. The juodel I is 
somehow related to the models studied inHmfca, and 
the model II is the Kondo lattice model considered inH 
with ferromagnetic coupling. Applying Abelian and non- 
Abelian bosonization methods, we investigate low-energy 
properties of these models, and find that they show quite 
different behaviors at the fixed points. We also study the 
effects of randomness on these fixed points, which may 
play a crucial role for HDHS because the hole-doping in- 
evitably induces randomness to the system. By using 
replica trick and renormalization group analysis, we find 
that the quantum disordered state inherent in Haldane 
systems suppresses the Anderson localization for weak 
randomness. 

We first consider the model I without randomness. The 
Hamiltonian is given by 

H = -t 4a c i+i<y + h - c - + U Y1 c \\ c i]4[ c il 

to 

-t dtydi+ier + h.c. + U dl^d^dl^dii 

+jj2 s c,i-Sd ti , (i) 
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where Cj jCr and cj CT etc., are the annihilation and creation 
operators of electrons, U is on-site coulomb interaction 



in each chain, and S c — d tr ai pCi p/2 with the Pauli 
matrix cr, etc. Here the last term is the Hund coupling 
interaction with J < 0. At half-filling, in the strong 
coupling limit J — > — oo, the model reduces to the S = 1 
Heisenberg chain which possesses the Haldane gap. Even 
for small \J\, the system has a finite gap, antLhelongs 
to the universality class of the Haldane stateoEM. We 
study low-energy properties of this model away from half- 
filling by applying standard Abelian bosonization meth- 
ods. We first linearize the dispersion of electrons and 
express the electron operators in terms of boson fields 
4> ca and Oca-, etc. which satisfy the_papmutation re- 
lation [cf> ctT (x), ca {x )] = iO(x — x )E3o. Introducing 
the boson fields, <f>± — (4> c ± (j>d)/y/2 for the spin sector 
and (jr± — {(j)^ ± 4*4 ; )/V2 for the charge sector where 

4>a = (<AaT - 4>al)/V2 and 4>{ c) = {4> a ~\ + 4> a i)/V2 with 
a = c,d, and their canonical conjugate momenta n± and 

(c) 

IL± , we can write down the bosonized Hamiltonian as 
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J dx(Js cos V^7r</>_ + 2 J4 cos V4~7r6L 



+ J5 cos V4~7n/> + ) cos y/Anc^ , 



(2) 



where v^f 1 , v v , K^f 1 and K v are Luttinger liquid param- 
eters in each sector. Here we dropped irrelevant terms 
with oscillating factors and irrelevant intra-chain back- 
ward scattering terms. We also omit Umklapp terms 
because we are concerned with the case away from half- 

(c) 

filling. The charge mode, cfr + , is completely decou- 
pled and hence it is described by U(l) Gaussian the- 
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ory with central charge c = 1 (Tomonaga-Luttinger 
liquid)!!.]. Initially in the renormalization procedure, 
Ji — J2 — J3 = J4 — J5 = J- These coupling con- 
stants are renormalizcd in different ways. From simple 
dimensional analysis, it is easily seen that Ji-term is al- 
ways irrelevant. In the case of strong correlation limit, 
U — ► 00, the scaling dimension of the field, cosV^tt^- 
is close to 1/2. For small J, the dimension of J^-teim is 
thus smaller than 2 and generates spin gap in (f> + mode 

(c) 

and charge gap in (f>_ mode. In this case the mass gap 
is also open either in <p_ mode or 9_ mode. The state 
with the mass gap in </>_ mode corresponds to the Ising 
Neel ordered state with the non- vanishing order parame- 
ter lim| l _ :)Hoo (-l) 4 ^(5fS'|) = (cos y7r</> + ) (cos y/n</>-), 
whereas that with the mass gap in mode to the 
quantum disordered state which is a weak-coupling ana- 
logue of the Haldane state characterized by the hid- 



den string 
(cos ypK<\>+ 



lim| i _ 3 -|^ 00 (5f expOX^i-i SfiSj) = 
The mechanism for the spin gap-for.- 
mation is similar to that for the coupled spin chainsclcil. 
In the present case, the spjn=|gap state with massive 8- 
mode realizes at half-fillingolld, and hence it remains mas- 

(c) 

sive away from half- filling, because (j)_ mode is already 
massive in this parameter region. Then the fixed point 
is the metallic state with the spin gap d la Haldane. A 
similar spin-gap state was obtained for U-J or Hubbard 
chains coupled via transverse hoppingoQ. At this fixed 
point, dominant correlations develop in the inter-chain 
singlet pairing and "4fcp" CDW. The order parameters 
are, respectively, given by S s(x) = (c La -(x)dR- cr (x) + 
c Ra (x)d L - a {x)), Ocdw(x) = (c LrT d^,d Ra >c R<J ). The 
correlation functions of these order parameters show 

algebraic decay; (Oss(%)Oss(ty) ~ x~ 1 t 2K + and 

(Ocdw(x)O C dw(0)) ~ x- K + )/2 . In our case K { + ] < 1 
and hence the fluctuation for "4fcp" CDW is more dom- 
inant than that for singlet superconductivity. 

We now discuss the effects of random impurities 
on this fixed point. The effects of randomness on 
Tomonaga-I^ttinger liquid were investigated by sev- 
eral authorsEro. Here we apply renormalization group 
methods used by Giamarchi and SchulzcJ. We introduce 
the random impurity potential, 

H imp = ^ / dx[^(x)(c\ L c a n, + d ] aL d aR ) + h.c], (3) 
with the Gaussian distribution, 



exp(-L»" 1 / C(x)£(x)dx). 



(4) 



Here we omit forward scatterings because they can be 
incorporated into the shift of the chemical potential and 
do not affect the fixed point properties. In order to deal 
with the quenched disorder we use a replica trick. We 
consider the case for weak randomness and incorporate 
only the first order contribution in D^. Then there is 



no coupling between different replica indices, which will 
be omitted below. Applying a standard renormalization 
group analysis^, we obtain the scaling equations up to 
the lowest order in J and D^, 
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where we did not display the equations irrelevant to 
the following discussions, and omitted the renormal- 
ization effects due to J2 term, because they do not 
change the intrinsic results. Since the quantum disor- 
dered spin-gap state realizes in the pure system, we ini- 



tially have 2 — 



> 0, 2 



< 0, and 



2 - l/K- - #i c) > 0. We can see from eqs.(§) and @ 

(el 

that the values of K+ and K_ are reduced by the ran- 

(c) 

domness, and hence 2 — K + — K_ ' > holds in the pro- 
cess of the renormalization. Then from eq. (|J) , J 5 scales 
to the strong-coupling regime and gaps open in <p + and 
<f>_ , resulting in K+,K^/ — > 0. Moreover eq.(||) indi- 
cates the fixed point value of J3 is D^/{2 — K- — K^), 
and hence the terms including J 3 in eqs.(||) and (|l0|) are 
negligible to the lowest order in D^. As a result, eqs.(j|), 
(Q), and ( |iT)| ) become 



d£> 5 
~dT 
dJ 4 
~dT 
dK_ 
dl 
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(12) 
(13) 
(14) 



We also note that the correction of K- by randomness 
can be neglected in the right-hand side of eq. ( |i~2| ) up to 
the first order in D^. So, is solely determined by 
cqs.(|l3|) and ( pT| ) without the term proportional to D^. 
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Fig.l. Plots of l/K-(l) vs D^(l) for some initial values 

of J4. The initial values of 1/K-, D^, and are taken 
as 0.7, 0.01, and 0.6. 

Solving these equations, we have K_ — > oo, because the 
spin gap in mode opens in the absence of randomness. 
Thus from eq. (p"2|) , scales to and the randomness 
becomes irrelevant. These arguments are confirmed by 
the numerical results for renormalization flows calculated 
by eqs.fQ ~(0) (Fig.l). Note that the suppression of 
is due to the spin gap formation by J4-term, that 
is, the presence of the non- vanishing disorder parameter, 
(cos-JttOJ) . Therefore we can say that the Tomonaga- 
Luttinger liquid state in the charge sector is protected 
from the transition to the Anderson localization by the 
quantum spin disorder. The physical reason of the sup- 
pression of the Anderson localization may be understood 
by observing that the presence of the Haldane gap re- 
duces the number of massless excitations and hence has 
a tendency to effectively suppress backward scatterings 
due to impurities. 

We have checked that for sufficiently large initial values 
of and small initial values of J4, scales to a value 
larger than unity before J4 scales to the strong-coupling 
regime, and our weak-coupling treatment for randomness 
breaks down(see also Fig.l). Although in this case we 
do not know the fixed point properties, it seems plau- 
sible to expect that sufficiently strong randomness may 
destroy the Tomonaga-Luttinger liquid state and bring 
about the transition into the Anderson localization state. 
It may also leads to the destruction of the quantum dis- 
ordered spin-gap state. Although our discussions here 
for localization-delocalization transition rely on weak- 
coupling renormalization group methods, we think that 
the qualitative features should not be changed even if we 
take into account higher-order corrections in J. 

We wish to note that the suppression of the An- 
derson localization for weak randomness is not due to 
the development of superconducting fluctuation as found 
beforec3t3, because in the present case the fluctuation of 
the "4fcp" CDW is more dominant than that of super- 
conductivity. The presence of the quantum spin disorder 
inherent in the Haldane state indeed prevents impurity 
potential from pinning the CDW. It may be a new type of 
depinning effect of the CDW. Hence, if a doped-Haldane 
system belongs to the class of model I, we can observe the 



power-law behavior in correlation functions more easily, 
owing to the stability of Tomonaga-Luttinger liquid state 
against randomness. 

We next discuss the model II. The Hamiltonian is given 
by 

H = -t^4* c i+i<r + h - c - + U Y1 C h c ^4i c n 

ia 

+ Jd ^ Sd,i ■ Sd,i+1 + J ^ Sc,i ■ Sd,i, (15) 

i i 

where Jd is the anti-ferromagnetic exchange interaction 
(Jd > 0) and other notations are the same as those of 
eq.(|l]). The charge degree of freedom for d-electrons 
is completely frozen in this model, which may approx- 
imately describe the system in which the energy levels 
of electrons composing the S — 1 state are largely sep- 
arated. In the strong-coupling limit at half-filling, this 
Hamiltonian also reduces to the S = 1 Heisenberg model. 
For small |J|, it can be a weak-coupling analogue of the 
Haldane system as in the case of model I. Note that for 
J > 0, this model coincides with the Kondo lattice model 
with the nearest-neighbc* exchange interaction between 
localized spins studied intj 

In order to properly describe symmetry of the sys- 
tem away frcoaiialf-filling at the fixed point, non-Abelian 
bosonizationOlIa is more suitable than Abelian bosoniza- 
tion because two chains decouple at the fixed point and 
SU(2) symmetry of each chain should be retained as we 
will see below. By applying non-Abelian bosonization 
formula to the spin sector, we can write down the Hamil- 
tonian in a continuum limit for the case away from half- 
filling as 



H = H C + H. 
H r = I dx 



v (c) i, {c) K {c) 

^y(^#) 2 + ^f^n^ 



2K, 



f 2irv c 

H s = J dx [J C L ■ J cL + J cR ■ J cR\ 



dx — —^[JdL ■ J d.L + J dR ■ JdR] 



-Jm J dx[J c l ■ J dL + J cR - J dR] 



+J ir / dx[J cL ■ J 



dR 



' cR 



< dL\ 



(16) 



where J c l and JdL, etc. are the current operators for 
spinons of c- and (i-electrons, which satisfy SU(2) Kac- 
Moody algebra, and initially we have J m — Ji r = J. 
Here we dropped irrelevant oscillating terms. In general, 
spinons of c- and d-electrons have different velocities, v c 
and Vd- From the operator product expansion of level- 1 
SU(2) Kac-Moody algebra, we obtain the scaling equa- 
tions for couplings J m and Ji r , 



dl 



= 0, 



dl 



72 



2tt(v c + v d ) ' 



(17) 
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Thus the last term of H s is marginally irrelevant for 
J < 0. The J m -term just renormalizes the velocities 
of spinons and does not break SU(2) symmetry in each 
chain. Hence at the low-energy fixed point there exist one 
massless mode in the charge sector described by c = 1 
Gaussian theory and two massless modes in the spin 
sector .described by level-1 SU(2) Wess-Zumino-Witten 
theorjO, and consequently there is no spin gap away 
from half-filling. Therefore doping holes in this model 
drastically changes characteristic properties of spin ex- 
citations, in contrast to the model I. We note that the 
existence of two massless spin modes is also contrasted 
to the Kondo lattice model with the antiferromagnetic 
coupling J > where only .one massless spin mode with 
a large Fermi surface existsc3u. 

At this fixed point, two kinds of spinons have different 
pseudo-Fermi surfaces. Therefore correlation functions 
involving spin excitations show singularities at two dif- 
ferent points reflecting the existence of two pseudo- Fermi 
surfaces. For example, the spin correlation function for 
the total spins, S(x) = S c (x) + Sd(x), is given by 



{S*{x)S*(0)) 



A a A ie l2kFX A 2 e inx 



+ c.c. (18) 



X" X" X 

where Aq,A±, and Ai are some constants, and the ex- 

(c) 

ponent of the second term is given by a = 1 + Kc ■ 
Here we omitted logarithmic corrections due to marginal 
operators. Thus the form factor of the spin correlation 
function may exhibit two structures in the momentum 
space at q = 2kp and n. 



The effects of randomness on the model II can also be 
investigated. In this model, two chains are decoupled at 
the fixed point. Then the problem is simply divided into 
two parts; the ID Hubbard model and the ID 5 = 1/2 
Heisenberg model with site randomness. For_thc former 
problem, according to Giamarchi and Schulza, the An- 
derson localization takes place for repulsive interaction. 
For the latter problem, if the on-site correlation between 
localized d-electrons is sufficiently strong, the site ran- 
domness does not affect the exchange interaction Jd so 
much, and hence the massless spin mode can survive, 
although the system may be of the glass state if we ad- 
ditionally include randomness for exchange interaction. 

Comparing the results for the models I and II, which 
exhibit quite different behaviors upon doping, we can say 
that characteristic properties of hole-doped Haldane sys- 
tems strongly depend on whether the energy levels of 
electrons consisting the 5 = 1 state are nearly degener- 
ate, or largely separated. p 

According to the experimental results of reffl, low- 
energy excited states below the spin gap appear upon 
doping, and localization occurs in transport properties. 
These seem to be consistent with the model II with ran- 
dom potential. We think that the model II may effec- 
tively describe low-energy properties of the doped Hal- 
dane system Y2_ a; Ca a ;BaNi05 for larger doping regions, 
though hole-doping has-heen realized in a way different 
from the present modelnu. 
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